Shell Model for Time-correlated Random Advection of Passive Scalars 
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We study a minimal shell model for the advection of a passive scalar by a Gaussian time correlated 
velocity field. The anomalous scaling properties of the white noise limit are studied analytically. 
The effect of the time correlations are investigated using perturbation theory around the white noise 
limit and non-perturbatively by numerical integration. The time correlation of the velocity field is 
seen to enhance the intermittency of the passive scalar. 
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I. INTRODUCTION 



The advection of a scalar observable 
linear PDE 



{x, t) by a velocity field v is described in classical hydrodynamics by the 
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V- V0 = nV^e + i. 



(1.1) 



If V is assumed to be solution of the Navier-Stokes equations in a turbulent regime and the Pcclct num ber Pe, which 
measures the ratio between the strength of the advective effects and the molecular diffusion k in (1.1), is large 
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{L and v are the characteristic length and advection velocity in the problem) and if a steady state is reached, an 
inertial range sets in where both the effects of the forcing / limited to the large scales and those of the molecular 
diffusion acting mainly on the small scales can be neglected. In the inertial range no typical scale is supposed to 
characterise the flow. As a consequence, the structure functions of the scalar field 



5p(r) = ([0(x- 



(1.2) 



display a power law behaviour in the inertial range with anomalous scaling exponents H{p) |l|. The word anomalous 
means that the exponents H{p) deviates from the linear behaviour predicted by a direct scaling analysis of (1.1). 

It was first realised by Kraichnan IQ] that anomalous scaling can be observed in the mathematically more tractable 
case of the advection by a random homogeneous and isotropic Gaussian velocity field, which is delta correlated (white 
noise) in time and has zero average and covariance in d dimensions given by: 
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The power law behaviour of the covariance mimics an infinite inertial range for the velocity field. The scaling exponent 
^wn is a free parameter characterising the degree of turbulence of the adverting field. The physically meaningful values 
range from zero to two. In the first limit the effect of the random advection is just to define an effective diffusion 
constant ||]. In the latter case the velocity increments are smooth as it is expected for a laminar flow. The choice 
i,wn equal to four thirds represents the scaling of the velocity field conjectured by Kolmogorov for the solution of the 
Navier-Stokes equation in turbulent regime. 

The hypothesis of delta correlation in time is of great mathematical advantage for it allows to write the equation 
of motions of the scalar correlations in a linear closed form. The evolution of each correlation in the inertial range 
is specified by a linear differential operator, the inertial operator, plus matching conditions at the boundary of the 
inertial range. The occurrence of anomalous scaling has been related to the existence of zero modes of the inertial 
operators dominating the scaling properties of higher order correlations ( y,|^,§,^ and for a recent review and more 



complete bibliography |Q] ) . The behaviour of the anomaly has also been numerically measured for the fourth order 
structure function versus the turbulence parameter S^^n H- However, to implement accurate numerical experiments 
still remains a difficult task. Therefore it turns out to be useful to use shell model as laboratories to test ideas and 
results related to the full PDE model (see [Q for a general introduction to the shell model concept). In [|l6|,|^ two 
different shell models advected by a delta correlated velocity field mimicking the Kraichnan model were constructed. 
Anomalous scaling was observed numerically and in the simpler case |17| it was proven analytically that the anomaly 
of the fourth order structure function is related to the anomalous scaling of the dominant zero mode of the incrtial 
operator. 

The passive scalar advection by a white noise velocity field is a useful mathematical model, but it still very far from 
a physical realistic velocity field which displays both time correlations and non Gaussian fluctuations. A first small 
step in this direction is made by investigating how the introduction of a time correlation in a Gaussian velocity field 
affects the statistical properties of the scalar field. 

In the present paper we introduce a time correlated velocity field in a shell model. This is done by substituting the 
white noise with the Ornstein-Uhlenbeck process which provides exponentially decaying time correlations (Section ||). 
We investigate the model both analytically and numerically. By means of stochastic variational calculus, which we 
shortly review in Appendices and 0, we show how to rewrite the equations of motion for the scalar correlations 
in integral non closed form. Such an operation allows the evaluation of the correction to the white noise inertial 
operator stemming from the time correlated velocity field. This procedure has the further advantage that it creates 
a non-ambiguous relation between the coupling terms for the scaling exponent ^„„ of white noise advection to the 
scaling exponent ^ of the equal time correlation of the time-correlated velocity field (Section p| ). 

The inertial operators can be expanded around the white noise limit in powers of an a-dimensional parameter which 
is interpreted as proportional to the ratio e between the time correlation and the turn-over time of the advecting field. 
We focus on the features of the steady state. There we assume that the averages over the Ornstein-Uhlenbeck process 
of all the observables are time translational invariant. As a consequence the inertial operators become linear up to 
any finite order in e. 

In the white noise case, when e is equal zero, we generalise the procedure first introduced in |l7| and we show 
that the scaling of the zero modes of the inertial operator of any order is captured by focusing on nearest-shell 



interactions. The equations are closed with a scaling Ansatz (section |V|) by postulating that the scalar field is "close" 
to a multiplicative process. Furthermore we perturb the closure scheme in order to extract the first order corrections 
in e to the anomalous exponents for different values of ^ ranging from zero to two. The prediction of perturbation 
theory is an e dependence (non-universality) of the exponents except for the second order H{2) (Section ^. The 
overall result is analogous to the one obtained in jlj] where a Gaussian time correlated velocity field is considered for 
the advection of the scalar field in ( |1 . iD : the introduction of time correlation is seen to enhance intermittency. The 
anomalies vanish smoothly in the laminar limit ^ equal two. 

To examine the validity of the results from the analytical calculations and explore the regime with long time- 
correlations (e 3> 0) we turned to numerical experiments. The occurrence of corrections to the anomalies predicted 
by the perturbation theory for small values of e was confirmed. However, strong non-perturbative effects sets in and 
are dominating when the expansion parameter becomes of the order of unity. 

II. THE MODEL 

The model is defined by the equations {m — 1,2, ..., N) 

[^ + Kkl ]0„(i) - <5i„/(t) = z[fc„+iC+i(t)<„(^) - A;™C-i(^X-i(0] (2.1) 

Um{t) = — ^/ dse^'^^-qmis) (2.2) 

f{t)^^ f dse-'^vis) (2.3) 
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where the star denotes complex conjugation, the rjm{tys and 'q{t) are independent white noises with zero mean value 
and correlation: 

('ym(i)<(s))=2(5„„(5(i-s) and (77(^)77* (s)) = 2J(t - s) (2.4) 

The boundary conditions are ^o = Qn+i = 0. The model might be regarded as a severe truncation of the equation of 
the passive scalar (|l.l|) in Fourier space. The field component 9m is the representative of all the Fourier modes in the 



shell with a wavenumber ranging between km = kg A™ and k^+i = ^o A'"^^. The parameter A is the ratio between 
two adjacent scales and it is usually taken equal to two in order to identify each shell with an octave of wave numbers. 
The energy transfer in a turbulent flow is conjectured to occur mainly through the interactions of eddies of the same 
size. As a consequence the interactions in Fourier space are assumed to be local. The "localness" conjecture [Q is the 
motivation for the restriction to nearest neighbours of the couplings among the shells. 

In the absence of external forcing and dissipation the total "energy" of the passive field is conserved: 

, J N 

d^^=d^El^™l' = f°^ /W^'^^O. (2.5) 
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Far from the infra-red and the ultra-violet boundaries (i.e. for 1 <C to <C iV) the conservation of energy is expected to 
hold approximately giving rise to an inertial range. Equations ( ^.2| ) and ( |2.3| ) describe the random evolution according 
to Ornstein-Uhlenbeck (0-U) processes respectively of the advecting and external force field. The 0-U process has 
differentiable realisations thus making the random differential equations with multiplicative noise which specify the 
dynamics of the scalar 9 independent of the discretisation prescription. 
The velocity correlations are for t > s 

{um{t)ul^{s)) = ^(e-^ - e-^). (2.6) 

In the limit of large t only the stationary part survive. The a-dimensional parameter e appearing in the definition of 



the 0-U processes (2.2) and ([2.3[) defines the strength of the time correlation in units of the typical times t,„. In the 



white noise limit one has 

hm {um{t)uUs)) = 2\vm?5C-^^). (2.7) 
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The factor two is reintroduced here with the proviso that the delta must be understood according to the midpoint 
prescription when its argument coincides with the upper boundary of integration as it occurs in practical computations 
in perturbation theory around the white noise limit. 

For any finite e ordinary differential calculus holds true: the consistency conditions yields a Stratonovich discreti- 



sation prescription when e is set to zero and the recovery of the white noise advection model of |17|. 

The information about the scaling of the correlations of the velocity field at equal times is stored in the constants 
Vm- We assume the power law behaviour 

_i 

\Vm\ oc fcm^ (2.8) 

Kolmogorov scaling is specified by ^ = 2/3 while ^ = 2 corresponds to a laminar regime. The t,„'s in (p]2) describe 
the typical correlation times for the random velocity field. A simple physical interpretation is to identify them with 
the turn-over times, i.e., with the typical time rates of variation through non linearity of the advection field on each 
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ocfc„'+^ (2.9) 



The scaling of the correlation times is then fully specified in terms of the parameter ^. It is worth to note that for 
any ^ less than two the r^'s are always decreasing functions of the wave number. 

The evolution of the scalar is determined in the inertial range by its complex conjugate. It is useful to introduce 
an unified notation for the 2A^ degrees of freedom. Calling Q = 9^0* and U = u^u* one has for the N shells: 

, 2N 2N 

-e„ = Y, [A^^p + Y, BlpU^] Q/3 + /^a,i + rSo^M+i (2.10) 

/3=1 7=1 

with 



m,/3 = ~Kkmdm,p 




N+m,P = ~l^k^dm,n 




^^13 = —ikm+l[Sfj,m+l Sa,N+m " 


~ <5/3,m 5a,N+m+l 


Q,/3™ = ikm+l[Sp,N+rn+l 5a,m ' 


- Sp^N+7n Sa,m] 



(2.11) 



where Latin and Greek indices range respectively from 1 to TV and from 1 to 27V. The set of matrices with constant 
entries B'* do not commut e with in each other and with the A matrix. The known sufficient condition (see for example 
[ jl5[ ) to have a solution of ( ^.10 ) in an analytic exponential form is therefore not satisfied. From the geometrical point 
of view non-commutativity means that the dynamics is confined on a manifold which turns into an hyper-sphere in 
C^ in the incrtial limit ( b.5|). 

The complex equations (|2.10| ) are invariant under phase transformations. Given two diagonal hermitian 2N x 2N 
matrices with time independent random entries 



T = diag(e 
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S = diag(e-'("*i+'^^\ ..., e-'('^"-i+<^"), 0, e'('^i+<*2), 
if is a realisation of the solution of the equations of motion then 

Te([/) = e(5[/) 



P«(0JV-1+0n) 



0) 



(2.12) 
(2.13) 

(2.14) 



is still a solution. The phase symmetry is the remnant of the translational invariance of the original hydrodynamical 
equations in real space |10[| . From the phase symmetry (2.14) it follows that at stationarity the only analytic non zero 
moments of the correlation are of the form 






(nr=ie™,eAr+™,) = (n^^iie^f) 



(2.15) 



In the inertial range such quantities display a power law behaviour. The diagonal sector of the moments whose scaling 
properties are specified by the exponents H{2uj) 



(j(2 



1,-H(2u. 



(2.16) 



is in the shell model context the analogue of the structure functions (1.2) of the original PDF model (1.1). The 
exponents if (2cj)'s are said to be normal if they can be derived from dimensio nal analysis. Under the assumption 
that a steady state is reached one matches the scaling of the inertial terms in (2.1) with a power law decay of the 
solution 



The resulting prediction is a linear behaviour of the exponents versus the order uj of the diagonal correlation: 



H{2cj)=u;{l-^] 



(2.18) 



The scaling argument (2.17) neglects completely the random fluctuations of the passive scalar field. Normal scaling 
holds if the statistics of the ^-field is Gaussian. Deviations from normal scaling are then related to the occurrence 
of intermittency corrections to the Gaussian statistics. A systematic account of the fluctuations is provided by the 
study of the equations of motion satisfled by the moments of the scalar field. 



III. EQUATIONS OF MOTION OF THE FIELD MOMENTS 



In the white noise limit, e equal zero, the Furutsu-Donsker-Novikov formula W and the delta correlation in time 
of the velocity ensure that the moments C'^^'^^ are specified by the solutions of closed linear systems p^,E^. In the 
presence of finite time correlations stochastic calculus of variations [[l8|jl9[ allows to write non closed integro-differential 
equations for the correlations. A typical functional integration by parts relation is: 



(F{eit))UN+rr.{t)) 



ds 



{UN+m{t)U,n{s))(^^}^^RaAt^s)Bl^e^{s) 



\ rfe„(i) 



(3.1) 



where Finstein convention holds for r epeated Greek indices. The matrix R is the fundamental solution of the homo- 
geneous system associated with ( 2.10 ). A heuristic proof of the stochastic integration by parts formula and of (O) is 
provided in Appendices ^ and y . 

Let us start with the second moment of the scalar field 



C^rn\t)^ {e,nit)QN+.-.it)) EE (0™(i)C(t)> • 



(3.2) 



^From the equations of motion ( ^.1C| ) one has 



,^^+2Kkl]C^^\t)-2^{{QN+m{t)f{t))}5m,l = 

+ 2'R{ikrn+l {UN+rn{t)QN+rn+l{t)QN+rn{t))} + 

-2k^^{l {UN+rn-lit)eN+m-l{t)eN+rn-l{t))}. (3.3) 



The integration by parts formula (3.1) gives 



[^ + 2Kkl]C^Ht) - 2S,n,in J ds (/(i)/(s)*) {RN+n..N+l{t, s)) 



"'m+l'''m / ds J^-'m (^: *) ~^ 
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~2 fc^,r„._i / ds (^»^^iW^^+--i(^)) jR^(^^)^(,^ ,) (3.4) 
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where to = 1, . . . , A^, 5R is the real part and: 

•^m (''J ^)~^N+m+l,N+m ; N+m,m+l\''^ ■^) ~ ^N+rn+l,N+m+l ■,N+m,mV'^ ^/ 
~^'^N+m,N+m;N+m+l,m+lV-T^) ~ ^ N +m,N+m+l\ N+m+l,m\'^T ^) (3.5) 

G!^+m.iv+n;W+P.,(^'*)=E(Q^+pW^^+"+i."(*'0)^o,Wn(5,0)e,(s)) (3.6) 

a=l 



When a steady state is reached the LHS of (3.4) can be neglected through the whole inertial range. The RHS specifies 
the inertial operator of the theory. A further simplification is attained in the limit of very large shell number. For 
any ^ less than two one has 

m|oC' Tj-fi m'\oo Tm 



independently of e. At equal times the resolvent matrix R reduces to the identity. /^From (2.8) and (|2.9|) it stems that 



Hence for m going to infinity the inertial operator is linearised in the form 



'^m+l""rn '''m — A . (3.9) 



RHS = 2- 



(Ci?|i - C^^') - 2 ^ {C^ - CI^A . (3.10) 



The slowest decay scaling solution compatible with a zero LHS is 

C^ oc r„j = fc„ , (3.11) 

In other words we have proven that the scaling of the second moment is normal as it coincides with the dimensional 
prediction (p.l^) . Moreover since the result does not depend on e it is universal versus the time correlation. It is 



worth stressing that the derivation of (3.11) requires that each of the terms appearing in (3.10) has separately a finite 
non zero limit for m going to infinity. The condition turns out to be not self-consistent when the same reasoning is 
applied to moments higher than the second. 

An important consequence of normal scaling of the Cm s is the Obukhov-Corrsin ||ll| , [l2|] law for the decay of the 
power spectrum r(fe) of the passive scalar if the Kolmogorov scaling is assumed for the advecting field: 



r(fc) = ^ E ((^«^:)') « k-W^^+\^y, = fc~i. (3.12) 



d 

A second interesting limit is when e tends to zero. Neglecting all non stationary contributions to the velocity 



correlations the RHS of (3.4) becomes 



Rub — 2kj^j^idm y^m+l ^m ] 2fc,„d„j_i I6„ L'm-l J + 

nt - rl 



+kld,n-ij dse 4;:i^sR^(^2^^(i,5). (3.13) 

If e is set exactly to zero the integral terms disappear and the white noise equations of [1^] are recovered. The 
information about the scaling of the velocity field is absorbed in the (i,„'s. In a pure white noise theory it is convenient 
to redefine the turbulence parameter as 

?»„ = 1 + |. (3.14) 

A Kolmogorov scaling of the vel ocity field corresponds to ^wn equal to four thirds which is also the value giving 
the Obukhov-Corrsin scaling in ( 3.12| ). The two definitions of the degree of turbulence coincide for ^ equal to 



two (Batchelor limit). It is natural to identify ^u,„ with the turbulence parameter of the Kraichnan model. The 
correspondence fixes the physical range of ^ between [—2, 2]. 

In the general case of the 2ci;-th even moment of the scalar C'-^'^-' one has: 

RffQ = V^ t{2lj;0) r<{2(^) I 

E2P d / p"^^^_-iI?F(2") (f s) + 

+ Y^ 2fc^^d,„,_i J^se-^^^ ^^^J^^^^l^^^_,_^Jt,s). (3.15) 

The multidimensional matrix /(^"'O) jg ^^le linear inertial operator of the white noise theory. The integrand functions 



■^m m -1 m (^i *) are given by the straightforward generalisation of (B.f:). The LHS, as above, is set to zero as 



rni ,...,mi 

far as the steady state features of the inertial range are concerned. Repeated integrations by parts in the large t limit 
generate a Laplace asymptotic expansion ||l^ of integral terms in the RHS, the coefficient of which are the derivatives 
with respect to s of the functions .F'-^"^ evaluated at equal times. When the steady state sets in we assume the latter 
quantities to be invariant under time translations for large t. Under such an assumption it will be proven in section M 
that the equal time derivatives are specified at equilibrium by linear combinations of the C'^'^^^^s. The effect of a small 
time correlation is therefore to generate new couplings of order e in the inertial operators. The observables we focus 
on are the scaling exponents. As discussed in the introduction, anomalies occur in the presence of non trivial scaling 
zero modes of the white noise inertial operators. It makes sense to relate the e dependence of the anomalous exponents 
to a perturbation of the scaling zero modes derived for e equal zero. A straightforward approach to the problem calls 
for the solution of A'^" linear equations. A further source of difhculty is that the exact determination of the zero 
eigenvectors of the inertial operators of any given order requires the matching of infra-red and ultra-violet boundary 
conditions. In the absence of an exact diagonalisation any analytical approach must rely on closure Ansatze to solve 
first the white noise problem and then to yield the corrections to the zero modes by linear perturbation theory. 

IV. THE WHITE NOISE CLOSURE 

In the present section we present a closure strategy to compute the i/(2u;)'s in the case of white noise advection. 
As shown in the previous section the second diagonal moment is normal and universal versus the time correlation. 
The first nontrivial zero mode problem is provided by the fourth order inertial operator /t^-'^' . In p^ it has been 
shown that the anomalous exponent p4 

H{A)^2H{2)-pi (4.1) 

can be extracted up to a very good accuracy from the solution of only two non-linear algebraic equations. The 
stationary equations for C'^' in the inertial range far from the infra-red and ultra-violet boundaries are given by 



7-(4;0) ^(4) 

2A2 = 

'rn '■rn — 1 'n '■/i — 1 'rn 'n 

J- ^(4) J i_r^(4) , 9^ ,Wi^n+l W^m-l 



^(4) ^(4) 

^va '^rn — 1 

One recognises two kinds of couplings in Im',n;p,q- 

• "Global" or "unconstrained" interactions: the indices p and q range respectively from to — 1 to to, + 1 and from 
n— 1 to n+ 1. The coupling are independent upon the relative values of m and n: in this sense they are referred 
as global. 



• 



"Purely local" int erac tions: they occur only for |m — n| < 1 and correspond to the terms proportional to the 



Kroenecker's S in (4.2) 



Anomalous scaling in the incrtial range is strictly related to the presence of such purely local interactions. Were these 
latter neglected the fourth order moment would have a normal scaling solution 

C^l^«^r^. (4.3) 

The idea is to capture the anomalous scaling by looking at the "renormalisation" of global couplings by pure short 
range ones. Disregarding the boundaries, the system is invariant under a simultaneous shift of the indices. Hence 
assuming a perfect index-shift invariance there are, for the m-th shell, only two independent equations where (5-like 
terms occur: 

= V /(4;o) c^i) 

/ J ra.m\p.q^^p^q 
p.q 



= V /(^'") , C(4' (4 4) 



/ J 7n,m—l;p,q p.q 
p.q 



The third equation involving a purely local interaction of the 77i-th shell with its nearest neighbours 

/ ^ m+l,m;p,q^^p,q' 



p,q 



is generat ed f rom the second of the (4.4)'s by a simple index shift. Therefore it is not regarded as independent. 
The pair (4.4) contains all the relevant information to extract the scaling of the fourth moment. It forms a closed 
system of equations independently on the shell number to, as one imposes scaling relations to hold within the set of 
"independent" moments of fourth order: 

fi(4) ^ 7-"/-'(4) (A r\ 

r^ _ ,j. l,--ff(2) W4) Ua-^ 

where the integer n is taken larger than zero. As in the analysis of the interactions the concept of independence stems 
from the assumption of index shift invariance: the moments of the form C^_„ „ are immediately reconstructed once 



(4.5) and (4.6) are given 






Let us analyse the closure Ansatz in more detail. The first equation (4.5) is a global scaling assumption of the 
"diagonal" sector of the fourth moment. Its justification lies in the very definition of an inertial range. The second 
scaling assumption relates the diagonal sector to the non diagonal one via a marginal scaling. It is the analogous in 
the present context of an operator product expansion (OPE) in statistical field theory p^. There renormalisation 
group (RG) techniques are able to describe the scaling behaviour of correlations of fields sampled at large real space 

7 



distances one from the other. If an observable requires the evaluation of a correlation including the products of one 
field in two points at short distance i.e. {4>(x — dx)(j)(x + dx)...) the RG procedure cannot be directly applied. The 
problem is overcome by an OPE or short distance expansion. The prescription is to rewrite the product via a Taylor 
expansion in terms of local composite operators sampled just in one point. Such a point is now well separated from 
all the other appearing in the correlation function. The original correlation is substituted by a set correlations such 
that RG applies provided an extra renormalisation, renormalisation of composite operators (RCO), is introduced. 
The latter is understood by observing that in our example the first term in the Taylor expansion gives 

4){x + dx)4){x — dx) ^ (f>{x) 

The mathematical meaning of a field is that one of an operator- valued distribution. The product of two distributions 
at equal points i.e. (t>{x)'^ requires a regularisation before the cut-off is removed in order to make itself sense as a 
distribution. This is the content of the RCO. Finally at leading order the relation between the renormalised quantities 
reads for the above example 

Mx + dx)cf>{x - dx)]« ...) ^ c{dx) { [cl^ixf] ^ ...) . (4.7) 

Roughly speaking the small real space separations are associated with the UV behaviour of the Fourier conjugated 
variable. In the shell model context the 9m are representative of the scalar field variation over one octave. The 
moments C^ m+n correspond to the Fourier transform of t he fo urth order structure functions of a homoge neou s and 



isotropic real space solution of the passive scalar equation (1.1). The meaning of the non-diagonal closure ( |4.q ) is to 
assume the long range (many shells) behaviour of the Fourier transform of the OPE coefficient c{dx) to scale inside 
the inertial range for large n independently on m. The constant x renormalises the value of the first shell where 
deviation from scaling occur. The analogy with the OPE is then summarised by writing 

0(x) = [e{x) - eml 

lim f d^xe^'^'^'^cidx) ^ xk:^"}^^ (4.8) 

\k\^oc J 



The insertion of the scaling Ansatz in (4.4) leaves a non-linear system in the unknown variables z and x. By applying 
the definition kn = A" one gets 



-1 - A-^(2) +2x{l + z A-«(2) ) ^ 

(1 + z) A-^(2) + a; ;2 (_i _ 3 A-»(2) A-2 ff(2) + z A-3 ff(2)) ^ 

which after a straightforward manipulation provides z as the physical root of a second order polynomial 



(4.9) 



_ 1 + 2 A-g(2) + 2 A-2g(2) ^ x-^"^^^ + Vl + 4 A-g(2) + 8 A-2g(2) - 6 A-^^l^) - 4 X-5H{2) + x-^"i^) 
^- 2(2A-2^(2)+A-3«(2)+A-4ff(2)) • ^ ' 

In terms of z the anomaly is 

P4 = 2il(2)-J^ (4.11) 

log A 

and it proves proves to be in fair agreement with th e va lues obtained from the numerical solution of the exact equations 



(4.2) ||17[] and from the numerical integration of (Pj) for all the values of the turbulent exponent ^ in the physical 
range [0, 2] (see also Fig. g). The sign of p4 is always positive: the effect of the anomaly is to decrease the diagonal 
scaling exponent. 

The procedure presented in detail for the computation of the fourth order exponent is straightforwardly extended 
to any higher order moment when one recognises that two crucial observation hold in general. 

• In the absence of pure short range couplings the normal scaling prediction holds true far from the boundaries 
for the zero modes of the inertial operators of any order 2uj. 

• For any fixed shell m there is a one to one correspondence between the number of independent equations and 
moments of order 2uj. 



In the case of C^^'^-' there are 2"^"^ equations: for any fixed reference shell mi the interaction with the second index 
7712 is affected by a pure short range coupling if the latter is equal or one unit different from mi i.e. only two 
possible choices and so on until the oj-th index is reached. On the other hand 2'^^^ is the number of exponents which 
characterise the scaling of the 2a;-th moment. The Ansatz is that the marginal scaling of the non-diagonal sector 
is fully specified in terms of the diagonal scaling exponents of order less than 2uj. By means of the "OPE's" one is 
able to close the zero mode equations in terms of 2'^^^ unknown rcnormalisation constants and H{2uj). The analogy 
with a field theoretical OPE goes to show that the need for an infinite set of constant does not necessarily imply the 
non-renormalisability of the real space theory mimicked by the shell model |2l| ] . 

More concretely the diagonal scaling exponent of the sixth moment [uj — 3) of the scalar field 



^m,n,p — (0W+mO„i9Ar+„6„6jv+pOp) = {\0m\ \9n\ \0p\ ) 

according to the above criterion requires four independent equations (see appendix O) 

y- 7-(6;0) ^(6) _ 



(4.12) 



r(6;0) 

-m.7n,'. 

r(6';o) 
r(6;0) 



q ^. ^m.7n^m\p^q^r^p.q.r — ^ 

y- ' ' A^-fl) 






.(6) 



q.r m,7n — l,m — l;p,g,r Pi9i^ 



p.q.r 'm,m-\-l^m — l;p,q,rPiQ 



(4.13) 



r(6;0) 



(6) _ 



The OPE inspired closure yields 






a 



(6) 

m-\-n ,m-\-n ,m 



•^1 '"n-l ^m,n 



•ni~{-7i~{-p,m-{-n,7n ^ p—1 n — 1 m,m,r. 

Me) ^ ,-H(2)^(6) 



(4.14) 



Inserting the "OPE" in ( 4.13| ) one gets into the algebraic system for the unknown rcnormalisation constants (xi, a;2, X3) 
and the diagonal scaling factor z. 







-1 + A" (-l + 3zxi) + 3a:3 = 

-X^" z xi + X'^'^+P" z^ xi ~ 2 z (a;i-2a;2) + A" (1 + z (-7x1+4x3)) 

A" (l + 4a;i-6a;3) + A2" (4za;2 - 2x3) - .X3 = 

X^^+P" ZXl+X'^Z {Xl ~ '3X2) - ZX2 + X^^+P" Z^ X2 - X^" Z {X2~X3)+X'^°' {-AZX2+X3) =0 



(4.15) 



After some algebra (4.15) reduces to a single third order polynomial specifying the physical root of z. It is worth to 
remark that from the functional dependence of the coefficient of (4.15) the exponent H{6) depends upon the anomaly 
of i?(4). Once again the anomaly evaluated from 



P6^'3H{2) 



logz 
log A 



(4.16) 



is in fair agreement with numerics (see Fig. |2|) for different values of f . 

In appendix ^ the same steps are performed in the case of the eight moment Cm,n,p,q- The analytical predictions 
for the anomalous exponents are summarised in Fig. nl In all cases the anomalies are decreasing functions of the 
turbulence parameter ^ vanishing smoothly when the laminar limit (^ equal two) is approached. The anomaly of 
the fourth order moment can be compared with the results of numerical experiments for the fourth order structure 
function of the Kraichnan model [g . There the adopted turbulence parameter is ^u,„ . For values of ^u,„ of order one, 
i.e., from the Kolmogorov scaling up to the Batchelor limit one indeed observes the same monotonically decreasing 
behaviour with values of the anomaly of the same order of those found in the shell model. For lower values of ^^„ 
the anomaly in the Kraichnan model display a maximum before decreasing to zero for ^u,„ equal to zero, i.e., when ^ 
tends to minus two. No sign of such behaviour is observed in shell model. The discrepancy might be an artifact of the 
shell model, which was originally designed to mimic the supposed local in scale character of the nonlinear interactions 
in a turbulent flow |lO[ , fails to describe a regime where strong non local effects become important. 

On a phenomenological level the energy transfer in the inertial range of turbulent field is related to the occurrence 
of a cascade mechanism as firstly conjectured by Richardson p2[ . The conservation of energy in the inertial range 



imposes that the forcing occurring on large real space scales is transferred to small scales (i.e. large wave numbers) 
before being dissipated. A mathematical description of the cascade is provided by multiplicative stochastic processes 
p3[ . Multiplicative modelling has been shown to account for most of the features observed in real and synthetic 
turbulence [ pipcll . In the present case the idea of a multiplicative structure is incorporated in the hypothesis that 
the scaling of the non diagonal sector of a given moment of order 2 uj is reconstructed once the scaling of the lower 
moments is known. Such an assumption together with the analysis of the couplings in the inertial operator of order 
2 Lu yields with fair accuracy the scaling exponents of the model without the resort to an exact diagonalisation of 
inertial operator. 

V. PERTURBATIVE ANALYSIS 

Let us now turn to the time correlated case. The idea is to evaluate the scaling behaviour of the dominant zero 
modes of the inertial operators ( ^.15| ) linearised up to first order in e by perturbing the white noise closure Ansatz. 

The first order corrections in e to the inertial operators are obtained by truncating the integration by parts to the 
terms linear in e 

9h'^ d i p ^^^_si?7:"(2'^) (f a\ — 

^"-mi + l^mi / *- 1 •"-'^mi,...,mi,...,m^\^^^) ~ 

^2eX^-^^^Ti^^-l^^^_^Jt,s)\s=t + Oie^ee-^). (5.1) 

The use of ( |3.9| ) in the RHS stresses that the effective a-dimensional expansion parameter is eX^: the range of 
reliability of first order perturbation theory is compressed to e < A~^/10. As mentioned in section III, in the limit 



t going to infinity, one expects two time quantities to be stationary. In such a case the derivative with respect to 
the variable s can be interchanged with the derivative with respect to t and one can use the equations of motion to 



evaluate (5.1). A direct differentiation with respect to s is consi sten tly t ake n with respect to the system of stochastic 



differential equations conjugated by time reversal of equations (|2.l| ) to (2_^). The latter operation in general requires 
the knowledge of the probability density of the forward problem. In the stationary limit the time reversal operation 
for the 0-U process reduces to the inversion of the sign of the drift term as in the deterministic case. After a slightly 
more lengthy algebra the result is equal to the differentiation with respect to t with opposite sign. 

The computations in the general case are very cumbersome (see appendices O O and 0). It is convenient to 

(2) 

exemplify the procedure in the simpler case of the second order correlation. There are four contributions to diJ-'m ■ 
d ^(2) 






..^N^^^l^N+m+l;N+'m,m\^T^)\t=s ~ '^^m+l'-'m \t) + \^ N +m\t)^ rait) ) 



'T;'^N+m,N+m;N+jyi+l,m + li^^^)\t=s = ^'^^m^m (*) + (^©Af+m+l (^©m+l (^) ^ ■ (5.2) 



By definition 



(eAr+,„(i)e,„(i)) = i|(|0„,(i)n 
ew+,„+i(i)e„,+i(t)) = i^ (|0™+i(i)p) 



(2') 

The time derivative of '^J-]^_i is derived by a simple index shift. As the terms non-diagonal in the resolvent R indices 
are zero, the second moment inertial operator is not by affected by first order perturbation theory. The result is 
not surprising. The second moment has only one free index. Hence at any order of perturbation theory only global 
coupling can be generated which are forced by the symmetries of the model to be consistent with a normal scaling 
of the zero mode. Moreover the 6m 's components of the scalar evolve only through the coupling with their complex 
conjugated Qiq+rnS'- their variation is a second order effect in e. The only possible non zero corrections are viscous 
and can be consistently neglected. 
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Let us turn to draw the general picture when lj is larger than one. Once again the phase symmetries ( ^.14 ) and 

the fact that when 'SiTm^...,mi,...,m,^ is known all other terms are yielded by index shift or exchange operations, 
prevent the corrections to the global couplings from affecting the scaling properties: the resulting "global" sectors of 
inertial operators have a normal scaling zero mode. This is in agreement with the observation made in p4[ |, where 
the dependence of the scaling exponents on the time correlation for generalised models of passive scalar advection is 
predicted to appear only through anomalies. The corrections to the purely short range couplings are therefore the 
relevant ones. They occur in two ways. On one hand new terms of order e show up in the purely self and nearest- 
neighbours interactions. On the other hand, terms proportional to Smi,m ±2 appear. The latter ones are the most 
dangerous for they in principle perturb the logic of the white noise closure by introducing new independent equations 
and henceforth the need for more renormalisation constants in the non diagonal sector of the moments. Nevertheless 
one can argue a priori in the spirit of the renormalisation group pq ], only the nearest neighbours interactions are 
relevant for scaling. Hence first order corrections can be obtained allowing an e dependence in the renormalisation 
constant of the white noise closure and determining the first order coefficient of their Taylor expansion. Moreover 
for uj larger than two such a strategy is already able to take into account the corrections due to the purely second 
neighbours interactions. 

Let us analyse in more detail the case of Cm'n- The white noise closure is perturbed by introducing an e dependence 
in the renormalisation constants 

C^m+n.m ^ ^i'') K-1 ^m.m- (5-4) 



The marginal scaling in the non-diagonal sector in (5^) is assumed to stay universal as it is for C*-^-* while the e 



dependence is stored in the pre-factor. The diagonal exponent is then determined up to first order as 

g(4,6) = ^+6A %^ f a^ (5.5) 

log(A) A2zlog(A) 

where z' is the derivative of ^ at e equal zero yielded by the perturbative solution of the system 

\ Z^p,q[^m,m-l;p.,q ^ ^ ^m,m-l-p,q "-^P.gV'j — U- 

The correction to pi due to time correlation increases the anomaly leading to a slower decay of the diagonal moment. 
For z' is negative (see Fig. ||) the overall anomaly is 






e\' 



A2zlog(A) 



(5.7) 



In the range of reliability of first order perturbation theory the effect is very small: for eA^ w 0(10^^) the prediction is 
a correction amounting to the three percent of the white noise exponent H{A). The perturbative scheme just proposed 
does not take into account the emergence of pure second neighbours interactions. In order to weight their relevance 



for the diagonal scaling and simultaneously to check the hypothesis of normal scaling for the marginal scaling in (5.4) 
one can relax the closure in order to encompass the equation 



/ Am.' 



m~2;p,q + ^ ^m,m-2;p,ijj ^p,g (^) " (5-8) 



which describe the independent, in the sense stated above, second neighbours interaction. Consistency with the white 
noise theory imposes the latter equation to decouple when e is set to zero. The requirement is satisfied if the closure 
is chosen in the form 

Cln\-n,m = a;(6)q(e) '"""""'' fc(„_i)-ff(2) C^^ (5.9) 

the pre-factor q{e) 2 enforces q{0) to be a function of the white noise renormalisation constants. Were the 

white noise closure exact it would fix the value of g(0) to one. 

In Fig. ^the q{0) is plotted versus ^: through all the physical range it stays close to one with a maximal deviation 
on the order of four percent for ^ equal to minus two. Moreover as shown in Fig. the time correlation induced 



correction to H(4,e) when (5.8) is included has the same qualitative behaviour and is quantitatively very close to 
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the nearest-neighbours prediction. The result is an a-posteriori check of the robustness of the closure approach. It 
confirms that first order corrections can be extracted within the logical scheme of the zero order one. It follows that 
the equations specifying the zero modes of the inertial operator acting on the sixth moment (see appendix Q can be 
closed by assuming: 

^^ m ~\~ n,rn~{- 71, m~\-n V^/ ^^'m,rii,m 

m~\-n ,rn~\' 71 ,7n J- V / ri— 1 ^^7n,m,m,m 

-7(6) ^ . x,-H(2),-H(4,e)^(6) 

^^7n-\-7i-\-p,7n-{-7i,7n *^o y^jf^p^i '^n— 1 ^^7n,m,?n,m 



The exponent iJ(4, e) is known perturbatively from (5.5) while H{2) is universal. With the same rationale (ap- 
pendix u^) one can evaluate H{8, e). 

In Fig. the analytic predictions for the corrections to the scaling exponents are summarised. In all cases the 
corrections are negative i.e. they carry a positive contribution to p2w- The corrections increase with w, the rate of 
the growth being slightly slower than the A(2cj) ex lv^lj — 1) A(4)/2 predicted in jlj] for time-correlation generalised 
PDE Kraichnan models. 

Within the range of first order perturbation theory the overall effect of time correlation is seen to enhance inter- 
mittency. An intuitive understanding of the phenomenon might be obtained by interpreting the time correlation as 
a mechanism to increase the probability of coherent fluctuations of the scalar field. The latter are rare events felt in 
the tail of the probability density of the scalar field as extreme deviations from the Gaussian behaviour of the typical 
events. 

VI. NUMERICAL EXPERIMENTS 

The resort to numerical experiments has a double motivation. On one hand they can be used to test the predictions 
from the first order perturbation theory. On the other hand they provide a broader scenario of the features of the 
model beyond the grasp of perturbative approaches. The first task is far from being easy because a quantitative check 
of perturbation theory requires measurements of the scaling exponents within an accuracy smaller than two percent. 

The main feature of the inertial range is the conservation of the scalar "energy" . From the analyti cal p oint of 



view this is seen in the non-commutativity of the terms associated with the multiplicative noise, B'^ in ( 2.10| ). This 



property rules out the use of a simple Euler scheme, which can be applied in the case of delta correlated noise. In the 



case of white noise advection the multiplicative structure of the noise (2.10) , which is interpreted in the Stratonovich 
sense, can be mapped into the corresponding Ito equations. The advantage is that the diagonal non zero average part 
of the noise is explicitly turned into an effective drift term [E6l . The non-diagonal terms in the Taylor expansion of 
the scalar field O are of the order three halves in dt, which are neglected in the Euler scheme. This procedure becomes 
meaningless for a time-correlated noise. There ordinary calculus holds and in the Taylor expansion of O both diagonal 
and non-diagonal products of the noise are of the same order in dt. 

Moreover the algorithm to be used must tend smoothly to a white noise limit, so that the same relative error is 
preserved independently on the value of e. 

Following Burrage & Burrage |E^ a reliable way out from the mentioned difficulties is to adop t the Trotter-Lie- 



Magnus formula to integrate the equations of motion to first order. For each time increment d t (|2.10 ) is solved in 
exponential form. Fast matrix exponentiation algorithms are provided by the package expokit [ p8[ . 

To generate the correlated noise, the exact method described by Miguel & Toral |29] is used. This method ensures 
that the noise is accurate down to the limit e ^ 0. 

The relevant time scale to measure the convergence of the solution is the slowest time scale in the system, namely 
the eddy-turn-over time of the first shell estimated as the maximum between eri and ri . As shown in Fig. |6| more than 
Nr = 10000 eddy-turn-over times are needed to achieve a converged solution for the sixth order structure function. 
The time step is set by the fastest time scale in the system, which is the one of the largest shell etm- The number of 
iterations needed to achieve convergence is then for e less than one: 



NtTM ^ -/Vr ^(M-l)(l-|) 



#(iterations) = -^^^^ = H^Lx^'^'-m^-i) (6.1) 

which shows that the number of iterations needed grows like 1/e, making it difficult to get close to the white noise 
limit using the same algorithm. 
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The scaling of the diagonal moments of higher order has been extracted by means of extended self-similarity [Q , 
where the p-th order structure function is plotted versus the second order one, which is assumed to be normal. The 
scaling is found as the average slope of the logarithmic derivatives in the inertial range. 

We considered a system with twenty- five shells with wave numbers increasing as power of A = 2, with viscosity 
K = 5 X 10~^. This choice ensures that there are several shells in the dissipative range. We focused on the results for 
£, equal two thirds (Kolmogorov scaling). 

In Fig. the "normalised" structure functions (16*^,1) km are shown. The quality of the scaling is demonstrated 
by the fact that the moments show scaling over a wide range of scales. 

A summary of the numerical experiments is given in Fig. g where the scaling exponents are plotted versus the order 
of the moments of the scalar field for different values of e. It is evident that the anomaly grows as the time correlation 
increases. 

When turning to the interpretation of the results in more detail, the uncertainty in the extraction of the scaling 
has to be kept in mind. For the sixth moment this uncertainty turned out to be on the order of a four percent. 
The changes in the scaling between different values of e is also on the order of a few percent. This seems to exclude 
a proper resolution in the numerics to compare the results with the analytical predictions from the perturbation 
analysis. However, the results for different values of e can still be compared with some confidence, as the relative 
uncertainty between the different runs is much smaller than the absolute uncertainty. This means that the slope of, 
i.e., the sixth order structure functions vs. e will be well resolved, while the absolute values can be shifted up and 
down a few percent. 

In Fig. the analytical prediction of the exponents is compared with the result of the numerics. The theoretical 
points are systematically below the numerical ones which is due to the absolute uncertainty as explained above. The 
slope is the same for the analytical calculation and the numerics, giving credibility to the results of the perturbation 
analysis. It should be noted that the effect of time correlation on the anomaly is quantitatively quite small even in 
the non perturbative range when e is equal to one (eA^ equal four). 

The global picture provided by the numerical experiments is that H{2u}) is seen to be a non linear function of e 
which, after rapid initial decrease in the perturbative range, displays a much slower rate of variation. An interesting 
question is whether there is a limiting value of the scaling of the structure function as e 3> 1 or not. However the 
quality of the numerics does not allow us to answer this question. 

VII. CONCLUSION 

We have presented a shell model for the advection of a passive scalar by a velocity field which is exponentially 
correlated in time. We developed a systematic procedure to calculate the exponents of the correlation of the diagonal 
moments (the structure functions). For the delta correlated velocity we find good agreement between analytical and 
numerical calculations up to the eight order. We presented an analytical perturbative theory to compute the correction 
to the scaling exponents due to the exponentially correlated velocity field. 

The occurrence of anomalies in the exponents of the diagonal moments of the scalar and their non universality versus 
the intensity eA^ of the time correlation, is related to the presence of pure short range couplings in the corresponding 
inertial operator which provide for non trivial scaling of the zero modes. In the absence of such short range couplings, 
as is the case for the second moment, normal scaling would take place independently on the value of eA^. 

The behaviour of the anomalous exponents in the non-perturbative regime was studied numerically. This was found 
to be a non linear monotonic function of eA^, decreasing with a rate much slower than in the perturbative regime. It 
is thus clear that the addition of the time correlation to the advecting velocity field enhances the anomalous scaling. 
The anomaly found in the present study is still much smaller that what is found when the passive scalar is driven by a 
turbulent velocity field driven by Navier-Stokes turbulence or by a shell model for the velocity field M. This indicates 
that the non-Gaussian nature of the real turbulent velocity field plays a significant role for the strong anomalous 
scaling observed for real passive scalars. 
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APPENDIX A: STOCHASTIC INTEGRATION BY PARTS FORMULA 

A heuristic proof of the stochastic integration by parts formula is provided. For a rigorous treatment see ||l^, [[l9[ . 
Let Ct be a stochastic process whose reahsations are defined as the solution of the Ito SDE 



xt = b{xt,t) + a{xt,t)r]t , Xt\t=o = xq 

where rjt is white noise. Let Q the stochastic process specified by 

it = b{xt,t) + eh{xt,t)a{xt,t) +a{xt,t)rit Xt\t=o = Xq 



(Al) 



(A2) 



For e equal the two SDE's coalesce: ([A2|) can be derived from (Al) under the variation of the white noise r]t 



rjt + h{xt,t). The integration by parts formula states that for any smooth functional / the identity holds: 



^JiQ)) J.=0 = (^fiCt)J^ dsh{Cs,s)^ 



(A3) 



where (...)^ denotes the expectation values with respect to the measure induced by Ct In order to prove it let us 
observe that the transition probability density for ([A2|) can be written formally as a path integral (Ito discretisation) : 



Pr,'{x,t\xo,0) 



S^{x,t\xQ,0) 



If one introduces the functional 



Vx e 



-Sa^,*ko,0)+/J dt'[^ii^^<^£/i(xj,,t')-4'i'K'.*')] 



dt' x't' - b{xt') . 
2 ^ a{x[,t') ^ 



M{Ct) = e 



lo dt' f\^l\f,''^ eh{cc',,t')-4-h\r',,t')] 



one has by construction 



|(M(CD/(CD)ci=o 



(A4) 



(A5) 



(A6) 



To each realisation of the solutions of ( A2) corresponds a mapping rjt —^ xt = x{t, ?7t, e). Hence the last equality can 
be rewritten as the white noise average: 



de 



(M(a:(i,r;t,e))/(a;(i,77t,e))),,^=0 



which implies (A3) when e is set to zero. The derivative 

d 



de 



fiOU^^^DQdcJiCt) 



(A7) 



(A8) 



is a Frechet derivative. The dynamics of the stochastic process DQ is linear once the realisations Xt of Q are known: 

(A9) 



yt = Dxt 

Vt =ytdxAb{xt,t) +(7{xt,t)'qt] +h{xt,t)a{xt,t) 



It is worth to note that for 6 = 0, a — h = 1 the integration by parts formula (A3) reduces to 

t{dwj{wt)) = {f{wt)wt) 



(AlO) 



which is the Gaussian integration by parts formula (see e.g. ]i|) applied to the Wiener process Af{0,t). 

The generalisation to a multidimensional complex case proceeds straightforwardly by introducing 2A'^ variational 
parameters {ci,£*}f^i and applying the definitions 



<Vm{t)Vnis) >^2SmnSit-s) 



(All) 



for the white noise correlations. 



14 



APPENDIX B: STOCHASTIC INTEGRATION BY PARTS FOR THE O-U PROCESS 

As in the above appendix we limit ourselves to the real case the generalisation to the complex case being straight- 
forward. Functional differentiation is formally derived from a Frechet derivative with h{xt,t) = S{t — s) where s is a 
parameter specifying the time when the white noise is perturbed. The variation is assumed to be non-anticipating 
(causal) : 



lim / ds'd(s - s') = 



(Bl) 



Let us consider the system of SDE's 



im = Ki{x) + ^ a;n,n{x)Cn{t) (B2) 



where c is the coloured noise: 



n=l 



t 

e 



Cn{t)= / ds' ^??„(s') (B3) 



Functional differentiation gives 

N N „t 



— {D'lXm) ^y^D^Xk [dkbmix) + dky^a^^nix) / ds' ^ "_ rjnis')] + ' (Tm.l{x) (B4) 

The functional derivative is fully specified when it is known its form at the time s when the variation of the white 
noise occurs. The latter is determined by the causality requirement 

^(Dtcnit)) ^ [dt^^ + ^Sit ~ s)]dn,i (B5) 

at e^/T„ e,/r„ 



which implies the variation of the coloured noise to be nonzero only immediately after the instantaneous kick 

t-s 

A'c„(t) = ^^Sn,i Vt > s (B6) 



By differentiating (B4) one finds 



, „ (DiXm) = ^ cTm.iix) S{t - s) + smooth terms (B7) 

dt^ Ey^i 



^From the last equation it stems that for t — s 



— {D^x„i)\t=s = —;=cr„i,i{x) (B8) 

dt e.^fTn 



Consistency with (B4) then requires that the variation of the a;'s associated with a non anticipating variation of the 



white noise at time s fulfils the initial condition: 

Dlx^{s)\t=s^^ (B9) 

The integration by parts formula ( [A3| ) for a smooth functional 0(x) 

{0{xt)cn{t)) ^ ds'^^y^(D''xA,0{xt)) (BIO) 

Jo <^V^ ^ \ / 



The variation is the solution of the linear problem (B4) of which we define R to be the fundamental solution. It 
follows 
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N N 



{0{xt)cn{t)) = ds ds' ^^^([a,, 0(xt)]i?,,„,(t,s')(T„^„(a:;)) 



1=1 m=l 

Finally inverting the order of integration one obtains 

-1+£L N N 



(0(xt)c„(i)) = / ds'- 
This proves the real version of formula (|3.l[) . 



2e 



J2 E ([^-' 0{xt)]Rl,„^{t,s')a„,,ni<)) 



1 = 1 m=l 



(Bll) 



(B12) 



APPENDIX C: THE FOURTH ORDER CORRELATION TO FIRST ORDER 



The inertial operator acting on the fourth moment Cm.n{t) is in the large time limit 



RHS ^ 4^:°!p,,CW - 2fc2^+id™ J ds e-^ jfTl^lit. s)- 



^kldra-lj dse ^^^-i-^T^^l^Jt,s) 

-kl+idn dse~'^—nj'i%{t,s) 



(4) 



+Kd„-i dse ".-1— 3fi^;^!^i_^(i,5 



ds 



(CI) 



The bi-dimensional matrix Irnln:p,q is the white noise linear inertial operator. The corrections to the white noise 
theory are generated by the time derivative at equal times of the integrand function SRJvi.m 

•^m,n(*: *)= (0Ar+m(i)BAr+„(t)e„(t)-RAr+m+i,Ar+,„(t, s)e,„+i(s)]) + 

— {Q N +m{i)Q N +n{i)Qn{t)RN +rn+l,N +m+l{t-, s)Qm{s) 
+ {QN+m+l{t)QN+n{i)Qn{t)RN+m,N+m{i-, s)Qm+l{s) 

— {'dN+m+l{t)QN+n{t)Qn{t)RN+m,N+m+l{t, s)Qm{s) 
+ {'dN+m{t)'dN+m+l{t)QN+n{t)Rn,N+m{t, s)0,„+i(s) 

— {'dN+m{t)'dN+m+l{t)'dN+n{t)Rn,N+m+l{t, s)Qm{s) 
+ {<dN+m{t)'dN+m+l{t)Qn{t)RN+n,N+m{t, s)Qm+i{s)\} + 

— {'dN+m{t)'dN+m+l{t)'^n{t)RN+n,N+m+l{t,s)Qjn{s)]) (C2) 

After a double integration by parts neglecting viscous contributions one gets into 
' ' 2 "^'P ^ 

P><3 
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A^ 



A^ A^ 



A^ 



7eA4 leX" TeA* 1 eX" 



H 
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The diagonal scaling exponent is derived up to first order is e resorting to linear perturbation theory. If pure second 
neighbours interactions are taken into account the constant q{0) is specified by 

l-(l + A-^^(^)+A-^(^))z(0) 
^(') = z(0) + z(0pA-3H(2) (C4) 

The result is approximately equal to one for all S, ranging between [0, 2]. The first order correction z'(0) is extracted 
from the solution of the linear system: 

(4 a2+«(2) + 4 A2 z(0)) x'{0) + 4 A2 x(0) z'(0) = 
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[A2 (-3 - A-«(2) _ A«(2)) ^(0) ^ ;^2-2H(2) ^(o)2j ^,^Q^ ^ 
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APPENDIX D: THE INERTIAL OPERATOR FOR THE SIXTH MOMENT OF THE CORRELATION UP 

TO FIRST ORDER 

Under the hypothesis that pure second neighbours interactions do not require new equations to specify the diagonal 
scahng for small values of e, there are only four equations describing how global coupling are "renormalised" by 
relevant pure short range interactions. Given the rri-th shell one has 
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APPENDIX E: THE INERTIAL OPERATOR FOR THE EIGHT MOMENT OF THE CORRELATION UP 

TO FIRST ORDER 



The set of independent equations is finally given as: 
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The closure is provided again assuming scaling for all the possible conditioned expectation values with respect to a 
given shell. It follows 
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FIG. 1. The prediction of the closure Ansatz for the anomalfes in the scahng exponents of the fourth, p4, the sixth, pg and 
the eight, pg, moments of the scalar field versus the turbulence parameter ^. In all cases the anomalies are decreasing function 
of ^ going smoothly to zero as the Batchelor limit ^ equal to two is approached. The anomaly of the eight moment is obtained 
as the numerical solution of a sixth order polynomial. 
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FIG. 2. The analytical prediction for the anomalous part of the scaling exponent compared with the result of the numerical 
experiments for different values of the turbulence degree parameter 5. A Kolmogorov scaling of the advection field corresponds 
to ^ = 2/3. The dash-dotted line represents the (dimensional) normal scaling prediction. The continuum line interpolates 
the exponents as obtained from the numerical experiment (squares). The circles are the analytical prediction from the closure 
Ansatz. 
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FIG. 3. the first order corrections to -H'(4) with (dashed hne) and without inclusion of second neighbours interactions are 
plotted versus the turbulence exponent. The inclusion of second neighbours couplings increases the intensity of the anomaly 
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FIG. 4. The renormalisation constant q{0) is plotted versus ^. It remains close to one through all physical range proving 
self consistent the conjecture of normal scaling for the non-diagonal sector of the fourth moment. The result stresses that the 
renormalisation of nearest neighbours interactions provide the an accurate framework to extract the scaling exponent. 
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-2 -1^0 1 2 

FIG. 5. The first order correction to -ff(4) (continuous line), H{6) and H{8) versus the turbulence parameter ^. In the last 
two cases the corresponding linear systems are solved numerically. In all cases the corrections are derived by perturbing the 
white noise closure "renormalisation" constants. The effect of time correlation is seen to add a negative correction to the scaling 
exponents highlighting an increase of intermittency. 
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FIG. 6. The convergence of the sixth order structure function for e = 1. Shown is (\d%S) (t) for m, = 8. The fast upwards 
changes and slow downward relaxations reveal the intermittent nature of the signal ^^=8 (^) 
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FIG. 7. An example of the scaling of the structure functions for e = 2.0. The plot shows the structure functions "normalised" 
by the fitted scaling km to make the scaling regime appear as horizontal lines. The lower line is for p — 1, and the upper line 
is p = 8. Each line is offset to make it possible to distinguish the lines from each other. 
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FIG. 8. The anomalous part of the structure functions H{p) — p/3 as a function of p for e — 0.01 to e = 2. Also shown 
are the points from the analytical calculation of the scaling for the white noise case. The lines correspond to (from the top): 
e — 0.01, 0.02, 0.10, 0.25, 1.0 and 2.0. The dashed line corresponds to normal scaling. 
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FIG. 9. The scaling of the 6.th order structure function vs. e\^ . The inset shows an enlargement of the perturbative range 
eX^ << 1 where the analytical prediction from first order perturbation theory can be compared with the numerical experiments. 
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